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基于 边界 元 的 超 定 法 计算 板 切 口 应 力 强 度 因 于 
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摘 要 : 板 切 口 尖端 的 奇异 应 力 场 影 响 结构 安全 ,需要 对 含 切 口 板 进行 强度 计算 。 文 章 首先 根据 中 
厚 板 的 弹性 控制 方程 ,利用 切口 特征 渐 近 展开 理论 ,推演 出 中 厚 板 切 口 的 奇 性 特征 方程 ,引入 插值 
抵 阵 法 进行 求解 ,获得 板 切口 的 奇异 性 指数 、 特 征 角 函数 及 其 导数 。 之 后 使 用 边界 元 法 ,求解 切口 
板 内 离 切口 尖端 等 距离 圆 缴 上 点 的 位 移 ,并 将 其 代入 切口 尖端 位 移 渐 近 展 开 式 ,形成 一 组 非 正定 方 
程 组 ,采用 超 定 法 求解 ,获得 渐 近 展开 式 的 幅 值 系数 , 据 此 进一步 计算 出 切口 板 的 应 力 强度 因子 。 
文章 分 别 计 算 了 板 切 口 在 对 称 荷载 \ 反 对 称 荷 载 以 及 复合 型 荷载 作用 下 的 应 力 强度 因子 。 
关键 词 : 板 切口 ;应 力 强度 因子 ; 超 定 法 ;边界 元 法 
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Evaluation of stress intensity factors of notched plate by an 
over-deterministic method basing on the boundary element method 


LIU Jingwen , CHENG Changzheng, HUANG Yifan, WANG Feiyang, LI Tengyue 


(School of Civil and Hydraulic Engineering , Hefei University of Technology, 230009 Hefei , China ) 


Abstract :The singular stress field near the vertex of a notched plate will affect the safety of structure. It is 
necessary to evaluate the strength of a notched plate. The singular characteristic equation of the notch in a 
plate is firstly established , basing on the elastic governing equation of a plate and the asymptotic expansion 
theory near the notch tip. The interpolating matrix method is then introduced to solve it for getting the sin- 
gularity orders, characteristic angle functions and their derivatives of the notch. The boundary element 
method is thereafter applied to determine the displacements of the inner points on the arc around the notch 
tip. The displacement values are substituted into the asymptotic expansion to form a set of non-positive 
definite equations , which are solved by the over-deterministic method to obtain the amplitude coefficients 
in the asymptotic expansion. The stress intensity factors are further given out basing on the amplitude coef- 
ficients. The stress intensity factors of the notched plate under the symmetrical load, anti-symmetric load 
and mixed-mode load are respectively investigated. 
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板 作 为 工程 结构 常见 形式 ,广泛 用 于 航空 、 航 
天 汽车 .土木 等 领域 。 针 对 板 的 力学 分 析 , 最 初 使 
用 的 是 忽略 了 板 横 向 前 切 作用 的 Kirchhoff 薄板 理 
论 上 ,认为 板 的 位 移 仅 有 挠 度 。 对 于 薄板 ,由 于 转角 
位 移 贡 献 较 小 ,该 理论 可 较 好 地 适用 。 当 板 厚 逐 渐 
增 大 时 ,横向 剪 切 作 用 不 可 忽略 , 板 的 转角 位 移 对 板 
的 变形 影响 较 大 ,需要 建立 新 的 理论 来 描述 板 的 变 
形 , 由 此 产生 了 Mindlin 中 厚 板 理论 '” 以 及 Reissner 
中 厚 板 理 论 ” ,这 二 者 的 不 同 在 于 针对 剪 切 力 的 计 
算 所 采用 的 系数 不 同 ,具体 数值 可 参见 文献 [2- 4] 。 
本 研究 使 用 Reissner 中 厚 板 理论 研究 板 弯 曲 问题 。 

板结 构 在 设计 和 使 用 过 程 中 容易 形成 切口 , 切 
口 尖 端 会 产生 奇异 应 力 ,导致 结构 强度 下 降 。 切 口 
奇异 性 研究 始 于 平面 问题 , Williams P! 应 用 特征 函数 
展开 法 给 出 了 平面 裂纹 尖端 附近 的 位 移 场 和 应 力 
场 沁 后 续 此 法 被 应 用 到 板 裂纹 问题 中 ,得 出 了 一 些 
AE T. Kirchhoff 薄板 理论 的 ,较为 典型 的 薄板 弯曲 断 
裂 间 题 的 应 力 强 度 因子 解析 解 "。 由 于 Kirchhoff 
薄板 理论 忽略 横向 剪 切 变形 ,不 能 获得 含 裂纹 板 弯 
果 彰 亚 型 应 力 强度 因子 ,无 法 建立 完整 的 板 弯曲 断 
Afi Dt 为 解决 这 一 问题 ,基于 Reissner 中 厚 板 理 
话 的 断裂 问题 研究 受到 了 人 们 的 关注 。 对 于 含 裂 纹 
无 限 大 板 弯曲 断裂 问题 ,学 者 们 逐步 给 出 了 裂纹 尖 
端的 应 力 分 布 ”1 ,以 及 某 些 特定 荷载 作用 下 的 应 力 
强 麻 因 子 结果 ” ,而 后 根据 Williams 特征 函数 渐 近 
展开 法 ,给 出 对 称 / 反 对 称 位 移 及 内 力 的 表达 式 "” 。 
对 证 含 切口 的 中 厚 板 的 断裂 问题 , 由 于 其 应 力 奇异 
性 哥 切 口角 度 相 关 , 相 较 于 裂纹 问题 ,其 求解 难度 较 
KO Qian 等 ”应 用 双重 圭 级 数 展开 的 方式 获得 了 
中 厚 板 切口 尖端 的 位 移 及 应 力 表 达 式 。Yao 等 
基于 辛 对 偶 求解 体系 同样 给 出 了 中 厚 板 切口 /裂纹 
的 位 移 及 应 力 表达 式 , 并 与 有 限 元 法 结合 获得 了 相 
应 的 应 力 强度 因子 。 

边界 元 法 “中 作为 一 种 高 效 的 数值 计算 方法 ， 
相 较 于 有 限 元 法 ,具有 降 维 的 特点 , 且 由 于 采用 积分 
的 形式 获得 待 求解 方程 ,不 会 导致 精度 损失 ,因此 精 
度 较 高 。 学 者 们 一 直 致 力 于 研究 如 何 应 用 边界 元 法 
求解 板 弯 曲 问题 。 首 先 应 用 边界 元 法 解决 了 薄板 弯 
曲 问题 251 ,随后 又 将 边界 元 法 应 用 到 中 厚 板 弯曲 
问题 ,与 Reissner 中 厚 板 理论 结合 ,产生 了 专门 求解 
中 厚 板 的 边界 元 法 :1 。 

应 用 Williams 特征 渐 近 展开 的 方法 ,将 切口 尖 
端 位 移 场 的 渐 近 展开 式 引 入 中 厚 板 控制 方程 ,可 推 
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演出 中 厚 板 切 口 奇 性 特征 方程 ”。 切 口 奇 性 特征 
方程 与 裂纹 奇 性 特征 方程 形式 一 致 ,是 一 组 超越 方 
程 ,然而 由 于 切口 角度 非 零 ,导致 难以 获得 解析 解 ， 
一 般 采 用 数值 方法 进行 求解 , Niu 等 ”提出 的 插值 
和 矩阵 法 可 以 有 效 求解 这 类 特征 方程 。 切 口 奇 性 特征 
研究 ,只 能 获得 反映 切口 奇异 性 强 弱 的 奇异 性 指数 
以 及 特征 角 函 数 。 运 用 奇 性 特征 分 析 结 采 , 计 算 切 
口 的 应 力 强度 因子 是 为 一 个 重要 的 人 研究 诬 题 。 目 
前 ,有 限 元 法 O 研究 切口 的 应 力 强 度 因 子 主要 
方式 为 :将 切口 尖端 物理 场 的 特征 函数 渐 近 展开 与 
有 限 元 结合 ,求解 渐 近 展开 幅 值 系数 ,最 后 得 出 应 力 
强度 因子 。 应 用 边界 元 法 求解 中 厚 板 应 力 强度 因子 
的 文献 主要 集中 在 裂纹 问题 上 '“” 。 但 应 用 边界 元 
法 分 析 仿 切口 板 的 问题 尚未 展开 。 

本 研究 基于 切口 尖端 的 特征 渐 近 展开 技术 ,将 
切口 尖端 的 位 移 渐 近 展开 式 引 入 中 厚 板 控制 方程 ， 
获得 板 切口 的 奇 性 特征 方程 。 利 用 插值 矩阵 法 进行 
求解 ,获得 板 切口 的 位 移 渐 近 展开 式 的 奇异 性 指数 、 
特征 角 函 数 及 其 一 阶 导 数 。 再 利用 边界 元 法 计算 切 
口 板 在 丛 载 作用 下 的 内 点 位 移 , 然 后 将 所 得 位 移 值 
代入 板 切口 尖端 的 特征 渐 近 展开 式 , 形 成 一 组 非 正 
定 方程 ,运用 基于 最 小 二 乘法 原理 建立 的 超 定 法 ” 
求解 ,获得 板 切口 特征 渐 近 展开 式 的 幅 值 系数 ,最 后 
根据 幅 值 系数 、 特 征 角 函数 及 其 一 阶 导数 计算 出 板 
切口 的 应 力 强 度 因子 。 


1 中 厚 板 切 口 奇异 特征 分 析 


考虑 如 图 1(a) 所 示 的 中 厚 板 切口 , 板 宽 W, K 
L, 板 厚 h。 切口 板 中 面 见 图 1(b) ,切口 两 径 向 边界 
分 别 位 于 9=0, 和 0=60, 上 ,切口 张 角 @=0, -0,。 
考虑 板 的 横向 前 切 变形 , 板 的 位 移 分 量 包括 两 个 转角 
METE wr, ab, 和 一 个 挠 度 位 移 分 量 w, 见 图 1(c) 。 

将 位 移 分 量 几 ar, w 代入 中 厚 板 的 弹性 力学 平衡 
方程 ,可 以 获得 用 位 移 表 示 的 中 厚 板 控制 方程 ”为 


o -v OW, ^y, 
d x: ru. CEW y) zo 
ðx 2 ðy 2 ðxðy D\ ox 
14v O04, 1-v ow, ow, C [ Ow 
; M b. | -w,] =0 
2  oxoy 2 Ox oy D\ðy i 
Fw Fw dp, i dp, - 
ov^ Oy or Oy 
(1) 


式 中 :v 为 泊 松 比 ;D 和 C 23 8) AA 255 t 13] JS 0093 
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iz 


用 H 


EO Od TIAE FAJT] 8539 3$ 


切 刚 度 , 且 
Eh | 5Eh 


"a20bea 7C 340 9) (2) 


其 中 五 为 板 的 弹性 模 量 


切口 张 角 
2—0,-0, 


(a) 板 的 几何 参数 
图 1 


(b) 切口 板 中 面 


(c) 位 移 示 意图 


中 厚 板 切口 


Fig.1 Medium thickness V-notched plate 


为 方便 进行 特征 渐 近 展开 分 析 , 将 直角 坐标 系 

及 的 中 厚 板 控制 方程 进行 坐标 变换 ,得 到 极 坐 标 系 
民 鸭 中 厚 板 控制 方程 , 即 

PW. ,1 Wr s 1-o OV, 1o Owe 

o" r or prp 2° o 2 r or00 


3-v OU, C ðw 
. í x =0 
2r | 00 D 全 i.) 
Pav l. CAN vj, lv dpo 


2 r r6 PE 90 P ar? 
ol ap 1 po 1-v Vu, 


2 r r prb o9 2 7 
C/1 w 
5 r 0 ^) x 
2 oy, | ð 
£» ， 1 EM LUE ra Ue ce : Us =0 
ðr r Or r 00 ðr r r 00 


(3) 
式 中 :水 wr 29 dr, b, 转换 到 极 坐标 系 下 的 分 量 , 指 
向 参见 图 1(c)。 
为 引入 边界 条 件 ,这 里 仅 列 出 中 厚 板 的 弯 矩 
M, M PIDH 0,, 它 们 与 位 移 之 间 的 关系 为 


ac e ZL + 


Q,Cr,0) se i Ta 


M,,(r,0) 一 一 


(4) 


考虑 自由 边界 条 件 , 即 切口 表面 无 外 力作 用 ， 
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可 以 表示 为 
M,(r,8) |, 5,5, 70 
M Cr 28) longa 0 (5) 
Q,(r,0) laa 20 


中 厚 板 切口 尖端 附近 的 位 移 可 以 表达 为 级 数 渐 
近 展 开 的 形式 LB 


y, = Yo 十 2 Arp, C AL) 


V, = Yo + Y Av S O, A) (6) 


N 


wW = wo 十 Y A,r™w(0,A,) 
KTI 


式 中 :yo Wo 和 w 是 极 坐标 系 下 切口 尖 点 的 位 移 分 
18:4, 为 幅 值 系数 ;A 为 奇异 性 指数 ,下 文 简写 为 
A sir (0,4) air CO A) AI 0 (0,4) REGE ftt ER, 


下 文 分 别 简 写 为 ,Ww 和 w;NN 为 级 数 渐 近 展开 式 截 
取 的 项 数 ,截取 的 级 数 项 数 越 多 ,位 移 的 计算 精度 越 
高 ,能 精确 表达 的 位 移 场 距离 也 越 远 。 将 式 (6) 代 
人 入 式 (3) , 令 r 的 同 阶 系 数 之 和 为 0, 可 以 获得 中 厚 
板 切 口 i 即 


ls Ly, : 
2 V. m e p- V, 十 Ay, z0 
NE EE cus IN zl-4 _ 
ELDER ; d, =0 
w +A°w=0 

(7) 
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式 中 :(…) 和 (…) 分 别 表示 对 9 求 一 阶 导 和 二 
阶 导 。 

将 式 (6) 引 入 式 (4) ,再 代入 到 边界 条 件 式 (5 ) ， 
可 以 获得 用 位 移 特 征 角 函数 及 其 导数 表示 的 边界 条 
件 , 即 


(ur, t, uU U,) Med -0 
(lr, — rs tA We) locas, =0 (8) 
w 人 -0 


应 用 Niu 等 ”提出 的 插值 矩阵 法 ,求解 边界 条 
件 式 (8) 下 的 中 厚 板 切口 奇 性 特征 方程 式 (7) ,可 以 
获得 切口 的 奇异 性 指数 、 特 征 角 函数 及 其 一 阶 导数 。 

表 1 给 出 了 中 厚 板 切 口 前 三 阶 应 力 奇 异性 指数 
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计算 结果 ,其 中 ,A1、A, 对 应 面 内 拉 伸 和 面 内 剪 切 所 
产生 的 工 型 和 开 型 应 力 奇 异性 指数 ,A; 对 应 面 外 剪 
切 所 产生 的 五 型 应 力 奇异 性 指数 。 奇 异性 指数 来 源 
于 特征 函数 渐 近 展开 ,可 参见 公式 (6) ,不 同 角度 的 
切口 奇异 性 指数 大 小 不 同 。 一 般 来 说 ,裂纹 的 一 阶 
奇异 性 指数 最 小 ,应 力 的 奇异 性 强度 与 奇异 性 指数 
大 小 呈 负 相关 , 即 奇异 性 指数 越 小 ,切口 /裂纹 尖端 
的 应 力 奇 异性 越 大 。 从 表 中 可 以 看 出 , 当 切 口 张 角 
a >180" 时 就 会 产生 工 型 和 亚 型 应 力 奇异 性 ,而 只 有 
切口 张 角 o > 240" 时 才 会 产生 开 型 应 力 奇 异性 。 在 
表 1 中 ,本 研究 结果 与 参考 解 ” 进行 了 对 比 ,发 现 
本 研究 结果 与 参考 解 的 吻合 度 能 达到 小 数 点 后 五 位 
有 效 数字 。 


表 1 ， 中 厚 板 切口 前 三 阶 应 力 奇异 性 指数 


Tab.1 The first three singularity orders for the V-notch in a medium thickness plate 
Al À2 À3 

aC) 
本 研究 方法 2g o 本 研究 方法 25 fg Uo 本 研究 方法 2g o 
360 0. 500 000 0. 500 000 0. 500 000 0. 500 000 0. 500 000 0. 500 000 
| 340 0. 500 427 0. 500 426 0.562 010 0. 562 007 0.529 412 0. 529 412 
320 0. 503 491 0. 503 491 0.638 186 0.638 183 0.562 500 0. 562 500 
300 0.512 222 0.512 221 0. 730 904 0. 730 901 0.600 000 0. 600 000 
280 0. 530 396 0. 530 396 0. 843 444 0. 843 440 0.642 857 0. 642 857 
260 0. 562 840 0. 562 840 0. 980 480 0. 980 475 0.692 308 0. 692 308 
240 0.615 732 0.615 731 1. 000 000 1. 000 000 0. 750 000 0. 750 000 
| 220 0. 697 166 0. 697 165 - - 0.818 182 0.818 182 
200 0. 818 696 0. 818 696 - - 0. 900 000 0. 900 000 
180 1. 000 000 1. 000 000 " 一 1.000 000 1. 000 000 


中 图 2 给 出 了 切口 张 角 a =300" 时 前 三 阶 奇异 性 
指数 对 应 的 特征 角 函 数 ,从 图 中 可 以 看 出 ,关于 vy,、 
bs 的 特征 角 函 数 与 天 于 w 的 特征 角 函 数 处 于 解 灶 


1.0 
-150 -100  -50 


0 M 
-150 -100 -50 0 50 100 150 
0/(°) 


(a) A, 对 应 的 特征 角 函 数 


状态 。 实 际 计算 时 ,w 与 ur, ap, 的 关系 仅 体现 在 切 
口 尖 点 的 2 个 转角 位 移 所 带 来 的 刚体 位 移 。 


特征 角 函 数 
Cy 
[1 


-0.2 
-0.4 


L4 


.0 
0 50 100 150 -150 -100  -50 0 50 100 150 


8/(°) 0/(°) 
(b) A, 对 应 的 特征 角 函 数 


(c) 4, 对 应 的 特征 角 函 数 


图 2 PERY HORE KZE (a -300?) 


Fig.2 Characteristic angle functions for the V-notch in a medium thickness plate (œ 2 300?) 
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2 中 厚 板 切 口 奇异 强度 计算 


对 由 式 (1) 控制 的 中 厚 板 弹性 弯曲 问题 ,可 以 
通过 Betti 功 的 互 等 定理 ” 建立 中 厚 板 的 边界 积分 
方程 , 即 


caup) + [Ti(p,9) ug9) dT 


= | 000D CU -1,2,3) (9) 
式 中 :p 和 9g 为 源 点 以 及 场 点 ;wi 为 位 移 分 量 和 面 
力 分 量 , 当 j=1,2 IE uj t; 分 别 表示 转角 位 移 分 量 
y, ab, 和 由 弯 矩 组 合 形成 的 面 力 分 量 , 当 j =3 时 ， 
u; Jt 分 别 表 示 挠 度 位 移 w 以 及 由 剪 力 组 合 形成 的 面 
力 分 量 ;cy 表示 几何 系数 ,一 般 不 需要 单独 求解 ,可 
以 通过 刚体 位 移 法 消除 | T, (pa) uo) 47. 积分 项 
噶 竟 奇异 部 分 ,获得 非 奇 异 部 分 与 几何 系数 之 和 ; 
Usz=7, 分 别 表示 位 移 基本 人 解 以 及 面 力 基本 解 ,具体 表 
述 式 可 参见 文献 [18] 。 在 求解 域 的 边界 划分 单元 ， 
获得 离散 后 的 边界 积分 方程 , 即 


Co p) ujCp) F È | TDA 


DÈ ODEDA (ij =1,2,3) (10) 


式 惠 :M 表示 划分 的 单元 总 数 ;w REO RRE k 
号 单元 上 的 位 移 以 及 面 力 。 遍 历 边界 节点 , 列 边 界 
RIO Jr f X (10) ,获得 系统 的 边界 积分 方程 组 , 引 
入 这 界 条 件 后 可 求解 出 边界 节点 的 未 知 量 。 

应 用 上 述 边界 元 法 ,计算 获得 板 切口 尖端 附近 
各 总 的 3 个 位 移 分 量 p, ape, 和 ,将 其 与 渐 近 展开 
式 结合 ,得 到 包含 切口 尖 点 位 移 popo wo 以 及 待 
求 幅 值 系数 4, 的 方程 组 , 即 


N 
UM = Yio 十 cosg 》 Ar^ V. (0,AÀ,) "S 
k=1 
N ~ 
sin V A,r^ y ,(0,A,) 
k=1 
N ~ 
W = o + sin Y Ar y, (0,A,) + 
k=1 
N ~ 
cosg 》 A,r" Pal O, A) 
k=1 


N 
w = -— porcos — Yorsing + wo + $ Arn w(0,AÀ,) 
kzl 


(11) 


第 39 卷 
HPEH ORMAN r 的 圆 弧 上 nn 个 内 点 的 位 移 
值 进行 求解 ,具体 取 点 方式 参见 图 3。 


选取 内 点 的 岁 


(a) 切口 板 边界 单元 划分 (b) 超 定 法 内 点 选择 
图 3 基于 边界 元 的 超 定 法 取 点 


Fig.3 Chosen points for over-deterministic method based 


on boundary element method 
在 每 一 个 选 点 上 列 式 (11), 形 成 系统 方程 
组 , 即 


b. E 9 45 d 
Hy 0 ! 0 fà E 
W, -rcosg -rsing, 1 fi e fi itn 
. . . . . wo 
. EE . . . . . . A, 
V, 1 0 0 fa Did fw " 
Vs 0 1 0 fa 2 fon P 
. 3 3 i 
W, -reosü, -rsing, 1 f, ‘e f. 
(12) 


其 中 :yw wvwi(k=1,2,…,n) 表 示 选 取 的 第 上 号 
点 的 三 个 位 移 分量 ;4;(i=1,2,…,NN) 表示 第 i 个 渐 
近 展 开 式 的 幅 值 系数 ;9 表示 第 上 号 点 与 x RE I8] 
的 夹 角 ;f% 表 示 第 号 点 在 位 移 j 分 量 上 的 第 i 个 渐 
近 展 开 项 。 将 式 (12) 简 写 为 如 下 的 矩阵 形式 , 即 

Sa3nxl = Czas (N43) L (N43) x1 (13) 
IP iesu ua RZN n MAAR IRL SC, vs ANAE 
展开 系数 矩阵 ;Tv 1 WERU ORAM ELS UT 
近 展 开 式 幅 值 系数 的 列 向 量 。 可 以 看 出 ,当选 取 的 
节点 数 n> (N+3)/3 时 ,方程 (13) 为 一 组 超 定 方 
程 。 常 规 求解 线性 代数 方程 组 的 方法 失效 ,这 里 采 
用 基于 最 小 二 乘法 的 超 定 法 求解 ,其 超 定 解 为 


T-(C'C)"C'g (14) 
应 力 与 内 力 的 关系 式 为 
oy = EM, sra = Maur m zs -4z )Q, (15) 


中 厚 板 切口 的 应 力 强度 因子 计算 式 ”为 
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HME, Ap. 
Ki(z) = V2mlimr “ao(r;0,z) 
Katz) V2mlimr "^r, (r,0,z) (16 


Ki(z) = V2mlimr "p rar) 
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求 出 幅 值 系数 后 ,将 式 (15) 代 入 式 (16), 可 以 
得 到 用 幅 值 系数 、 特 征 角 函数 及 其 导数 表示 的 应 力 
强度 因子 , 即 


Ki (2) = -A, PEDE, (ON) d, 0,2.) 0 400,2) 
KyG) = =A, ED =v) (4 (0,:) -Bal0, As) + hs Wal0,A2)] (17) 
Ka (2) 2A, STC? Lag) 0,4.) 


从 式 (17) 中 可 以 看 出 ,Ki (2) 和 Ki (z) 的 最 大 


值 出 现在 板 的 上 下 表面 , 即 z = +h/2 Rb, Ky Cz) B] 


大 值 出 现在 板 的 中 面 , 即 z =0 处 。 对 中 厚 板 切 口 而 
言 ,3 个 应 力 强 度 因子 的 最 大 值 应 重点 关注 ,分 别 为 


K; = m Gp (0,4,) +y,(0,A1) + vA y,(0,A1)] 
Ki = -42 Sp - v) [p (0,A;) -yp (0,45) *À; i4 (0,A,)] (18) 
Ky A, STC i (0,A;) 


5 将 基于 边界 元 的 超 定 法 计算 得 到 的 幅 值 系数 ， 
导入 奇异 特征 分 析 得 到 的 奇异 性 指数 ,特征 角 函 数 
态 基 导数 代入 式 (18) ,可 以 计算 出 中 厚 板 切口 的 最 
大 应 力 强 度 因子 。 


3$ 例 


本 研究 依据 中 厚 板 切口 /裂纹 位 移 级 数 渐 近 展 
开 起 ,以 及 边界 元 法 计算 所 得 的 奇异 点 附近 指定 半 
做 正点 的 位 移 值 ,构成 待 求解 方程 组 ,并 使 用 超 定 法 
求解 超 定 方程 ,获得 幅 值 系数 ,进而 求解 应 力 强度 因 
子 。 首 先 给 出 含 中 心 裂 纹 板 在 对 称 弯 矩 荷载 作用 下 
的 应 力 强度 因子 结果 ,并 与 已 知 参考 解 进 行 对 比 。 
之 后 将 本 研究 方法 应 用 在 含 切口 板 受 各 类 荷载 作用 
下 的 应 力 强度 因子 与 取 点 半径 的 关系 研究 。 

1) 对称 弯曲 含 中 心 裂纹 方 板 ( 例 1) 


| 
| 
| Tq 
Lo 


M 
图 4 含 中心 裂 纹 板 受 弯曲 荷载 
Fig.4 A plate with a central crack subjected 
to the bending moment 
R2 ”弯曲 板 中 心 裂纹 应 力 强度 因子 Ki 
Tab.2 Stress intensity factors K, of a central crack 


in a plate under the bending moment 


图 4 所 示 的 方形 裂纹 板 AK L 2100 ,裂纹 长 度 
21 = 10 , 弯 矩 M=1, 四 边 简 支 。 弹 性 模 量 E =2 x10 ， 
泊 松 比 »=0.3。 计 算 不 同 板 厚 下 裂纹 右 尖 站 应 力 强 
度 因子 ,并 与 文献 [11] 和 文献 L26] 的 参考 解 进行 对 
比 , 结 有 果 列 于 表 2。 

从 表 2 可 以 看 出 ,本 研究 计算 得 到 的 应 力 强度 
因子 与 参考 解 的 误差 小 于 1% ,说 明了 本 研究 方法 
的 有 效 性 。 


h/l 本 研究 方法 Sg n 参考 解 二 [2 
1 0.7393 0. 746 2 0.744 6 
1.5 0.784 1 0.783 0 0.7810 

2 0. 820 4 0.8172 0.8150 


2 对 称 弯曲 切口 板 ( 例 2) 

图 5(a) 给 出 了 对 称 索 矩 作用 下 切口 板 受 力 示 
ARA] , 板 宽 歼 =50, 板 长 =100 ,切口 深度 为 /, 令 1= 
W/10, AB 边 与 CD 边 设 为 简 文 ,其 余 边 均 为 自由 。 
z 轴 根 据 右 手 定 则 确定 ,图 中 没有 画 出 。 为 获得 归 一 
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化 应 力 强度 因子 值 ,将 弯 矩 M 设置 为 1, 弹性 模 量 厂 
也 设 为 1, 泊 松 比 v 设 为 0.3。 


F 4. 
L 
3g 
E 


y 
Q 
ae 


t. 


(a) 对 称 弯 


(b) BATER Æ ih 


图 5 切口 板 受 弯曲 荷载 
Fig.5 A V-notched plate subjected to the bending moment 


四 根据 网 5(a) 所 示 的 加 载 模式 ,切口 板 在 受 载 后 


0.40 


—e— h=8.75 
—-— h-11.25 


—— h-1.5 
—— h-10 
一 一 h-12.5 


0.10 
ol 02 03 04 0.5 06 07 0:8 09 1O 11 


" 
(a) a-270* 


213 4530 35 
将 产生 沿 Oxz 平面 对 称 的 变形 ,因此 只 存在 工 型 应 
力 奇 异性 ,切口 的 应 力 强度 因子 也 只 有 Ki。 图 6 给 
出 了 利用 基于 边界 元 的 超 定 法 计算 得 到 的 应 力 强度 
因子 K 计算 结果 随 取 点 半径 > 的 变化 规律 。 

从 图 6 中 可 以 看 出 , 当 a 较 小 时 (如 w 270?) , 
取 点 需要 选择 相对 较 大 的 半径 才能 获得 稳定 的 应 力 
强度 因子 计算 结果 。 当 a 较 大 时 (如 œ =300°) , 计 
算 所 得 的 应 力 强度 因子 受 取 点 半径 的 影响 小 。 这 是 
由 于 当 a 较 小 时 ,用 于 计算 Ki 的 一 阶 奇 异性 指数 较 
大 ,一 阶 主 导 项 与 剩余 高 阶 项 之 间 的 数值 差距 较 小 ， 
对 主导 项 的 计算 精度 影响 较 大 。 从 图 6 还 可 以 发 
现 , 应 力 强 度 因子 计算 结果 的 稳定 性 与 板 的 厚度 之 
间 存 在 比例 关系 , 即 板 的 厚度 越 厚 ,计算 所 得 应 力 强 
度 因子 的 稳定 性 越 好 ,这 印证 了 本 研究 是 基于 中 厚 
板 理论 进行 的 奇异 强度 分 析 。 当 板 的 厚度 较 小 时 ， 
应 当 换 用 薄板 理论 来 建立 薄板 奇异 特征 方程 进行 求 
解 。 此 外 , 随 着 板 的 厚度 增加 ,应 力 强 度 因 子 的 值 会 
减 小 ,说 明 当 板 越 厚 时 , 越 不 容易 发 生 破 坏 。 


0.35 
—— h=7.5  —9—h-8.75 
0.30 —4— h=10 —— h-11.25 
一 一 h-12.5 


"Ol 02 05 04 05 05 07 UR 09 10 Ll 
n" 
(b) a=300° 


图 6 对 称 弯 曲 切口 板 的 K 与 取 点 半径 7 的 关系 


Fig.6 Variation of K, with r in the V-notched plate under symmetric bending moment 


3) BOSERKES H EL C9] 3) 

UEI 5 (b) Brzs B OSHPRSSABTEJ FARU O E 
型 , 板 将 沿 Oxz 平面 发 生 反对 称 变形 ,从 而 产生 [型 和 
亚 型 应 力 奇 异性 ,应 力 强 度 因 于 同时 包含 KM Km。 

图 7 给 出 了 基于 边界 元 的 超 定 法 ,计算 得 到 的 
反对 称 榴 曲 切口 板 的 下 型 应 力 强 度 因 于 天 gr 。 从 图 
中 可 以 看 出 , 随 着 切口 张 角 a 逐渐 减 小 ,应 力 强 度 因 
子 计 算 结 果 的 稳定 性 在 逐渐 减弱 。 随 着 板 厚 的 增 
加 ,应 力 强 度 因子 计算 结果 的 稳定 性 在 逐渐 增强 。 
板 的 厚度 越 大 ,其 对 应 切口 的 应 力 强度 因子 越 小 ,说 
明 破 坏 板 所 知 的 荷载 越 大 , 即 板 越 难以 被 破坏 。 

图 8 给 出 了 反对 称 芍 曲 切口 板 的 亚 型 应 力 强 度 


因子 Km 与 取 点 半径 了 的 关系 。 从 图 中 可 以 看 出 ,Km 
在 任意 半径 下 的 结果 都 趋 于 稳定 ,说 明石 型 应 力 强 
度 因子 的 计算 结果 不 受 取 点 半径 尺寸 影响 。 比 较 图 
8(a) 和 图 8(b) 会 发 现 ,在 不 同 切 口角 度 下 Kg 的 计 
算 结果 都 很 稳定 ,说 明 在 反对 称 加 载 模式 下 , 焉 型 应 
力 强度 因子 计算 结果 的 稳定 性 与 角度 无 关 。 因 为 无 
论 在 对 称 加 载 或 反对 称 加 载 模式 下 , 超 定 法 所 使 用 
的 边界 元 法 计算 得 到 的 位 移 值 中 , 挠 度 位 移 分 量 w 
都 占 主导 地 位 ,两 个 转角 位 移 分 量 办 和 办 相对 w 
较 小 。 因 此 , 挠 度 w 以 及 与 挠 度 w 相关 的 渐 近 展开 
项 的 计算 精度 相对 较 高 ,使 得 计算 所 得 的 焉 型 应 力 
强度 因子 非常 稳定 。 
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图 7 反对 称 弯曲 切口 板 的 Ki 与 取 点 半径 7 的 关系 


Fig.7 Variation of K, with r in the V-notched plate under anti-symmetric bending moment 
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各 4 ) 受 面 外 剪 切 荷载 的 板 切口 ( 例 4) 
二 图 9 给 出 的 受 面 外 剪 切 荷载 的 切口 板 , 其 几何 
尺 畦 以 及 材料 参数 与 例 1 相同 。 


PISS O ORR 
E 


图 9 受 面 外 剪 切 荷载 的 切口 板 
Fig.9 A V-notched plate subjected to the transverse shear 
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图 8 反对 称 弯曲 切口 板 的 Ky 与 取 点 半径 7 的 关系 


Fig.8 Variation of Ky with r in the V-notched plate under anti-symmetric bending moment 


AB 边 固 支 ,在 CD X156 99 UI fr r= 1, ER 
边 自 由 。 这 是 一 种 复合 加 载 模式 ,切口 尖端 将 同时 
产生 工 工 ` 亚 三 种 类 型 的 应 力 奇 异性 ,使 用 基于 边 
界 元 的 超 定 法 计算 ,可 以 同时 获得 3 个 应 力 强 度 
因子 。 

图 10 给 出 了 切口 张 角 a =300° 时 , 受 面 外 剪 切 
力作 用 的 板 切 口 的 应 力 强 度 因子 计算 结果 。 从 网 可 
以 看 出 , 受 面 外 剪 切 荷载 作用 切口 板 的 了 型 重型 
应 力 强度 因子 ,在 选取 不 同 半径 点 进行 计算 时 , 结 
保持 稳定 , 仅 有 开 型 应 力 强 度 因子 受 取 点 半径 影响 
较 大 ,但 随 着 板 厚 加 大 ,这 种 影响 也 逐渐 消失 , 即 当 
板 厚 较 大 时 ,Ki 也 不 受 选 点 半径 影响 。 因 此 ,在 复 
合 加 载 模式 下 ,基于 边界 元 的 超 定 法 可 以 准确 地 计 
算出 与 纯 工 型 或 正 、 焉 复合 型 切口 板 同样 稳定 的 应 
力 强 度 因子 。 
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(c) Ki 计算 结果 


图 10 受 面 外 剪 切 荷载 切口 板 应 力 强度 因子 (a — 300?) 
Fig. 10 Stress intensity factors of a V-notch in a plate 


under transverse shear ( & 2 300?) 


本 研究 基于 中 厚 板 的 控制 方程 ,采用 切口 尖端 
特征 渐 近 展开 理论 ,建立 了 中 厚 板 的 奇异 性 特征 方 
程 ,获取 了 中 厚 板 切 口 的 奇异 性 指数 以 及 对 应 的 特 
征 角 函数 。 再 应 用 边界 元 法 计算 切口 尖端 附近 内 点 
的 位 移 。 最 后 基于 奇 性 特征 分 析 结 果 以 及 内 点 位 移 
计算 值 ,形成 非 正 定 方程 组 ,采用 超 定 法 ,求解 出 渐 
近 展 开 式 的 幅 值 系数 ,进而 获得 中 厚 板 切 口 的 应 力 
强度 因子 。 

研究 发 现 中 厚 切 口 板 的 转角 位 移 特征 角 函 数 与 
挠 度 特征 角 枯 数 处 于 解 耦 状态 ,转角 位 移 对 挠 度 位 
移 的 影响 仅 体现 在 刚体 位 移 上 。 计 算 发 现 ,中 厚 切 
口 板 的 厚度 越 厚 ,应 力 强度 因子 计算 结果 受 边 界 元 
法 取 点 半径 影响 起 小 。 对 于 对 称 弯曲 切口 板 , 仅 存 
在 工 型 应 力 强 度 因 子 , 随 着 切口 张 角 的 增 大 ,其 计算 
结果 受 边 界 元 法 取 点 半径 的 影响 越 小 。 对 于 反对 称 
弯曲 切口 板 ,存在 开 、 亚 型 应 力 强度 因子 ,其 中 开 型 
应 力 强度 因子 计算 结果 的 稳定 性 , 随 着 切口 张 角 的 
增 大 受 取 点 半径 的 影响 越 小 ,而 亚 型 应 力 强 度 因子 
计算 结果 的 稳定 性 不 受 取 点 半径 影响 。 对 于 受 面 外 
荷载 的 切口 板 ,同时 存在 II 工 、 亚 型 应 力 强 度 因子 ， 
其 中 工 、 焉 型 应 力 强度 因子 ,在 选取 不 同 半径 点 计算 
时 ,结果 保持 稳定 , 卫 型 应 力 强度 因子 受 取 点 半径 影 
响 大 , 取 点 半径 越 大 ,计算 结果 越 稳定 。 

本 研究 所 提 的 边界 元 法 结合 特征 渐 近 展开 求解 
含 切 口 板 弯曲 问题 的 方法 ,可 以 有 效 计 算 工 程 中 含 
切口 /裂纹 板 受 弯曲 荷载 时 的 应 力 分 布 ,同时 也 可 以 
获得 切口 的 应 力 强度 因子 ,与 已 有 断裂 准则 相 结 合 ， 
可 判断 切口 板结 构 的 安全 状态 。 
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